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Abstract 

A one-way wave equation is a partial differential equation which, in some approxi- 
mate sense, behaves like the wave equation in one direction but permits no propagation 
in the opposite one The construction of such equations can be reduced to the approxi- 
mation of vl-s 2 on [-1,1] by a rational function r (s) = p m (s)/q n (s) This paper 
characterizes those rational functions r for which the corresponding one-way wave 
equation is well-posed, both as a partial differential equation and as an absorbing boun- 
dary condition for the wave equation We find that if r (s ) interpolates \/l - s 2 at 
sufficiently many points m (-1, 1), then well-posedness is assured It follows that absorb- 
ing boundary conditions based on Pade approximation are well-posed if and only if 
[m , n) lies m one of two distinct diagonals in the Pade table, the two proposed by 
Engquist and Majda Analogous results also hold for one-way wave equations derived 
from Chebyshev or least-squares approximation 
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1. Introduction 

The wave equation 

«« = + u yy (1 1 ) 

admits plane wave solutions 

«(*, y, t) = e’ (r<+ ^ +w) (12) 

for any that satisfy the dispersion relation 

r 2 = € 2 + V 2 (1 3) 

In the (£, rj) plane, the curves of constant r determined by this equation are concentric 
circles The phase and group velocities of (1 2) are identical and equal to 

v = (v x , v y ) = -- 2-) = (cos0, sm0), 9 £ [0,2tt] (1 4) 

r t 

Thus for each frequency r, (1 1) admits plane waves traveling with speed 1 in all direc- 
tions 

In certain applications, what is ideally wanted is a one-way wave equation or parax- 
ial equation that admits only half of these solutions in this paper, all those with v x < 0, 
l e 9 (E [y, 4p] This idea has been introduced over the years in various fields, and some 

early references are given in the Appendix to [Ta77] In the past decade and a half, 
one-way wave equations have come into heavy use in three areas First, they were intro- 
duced by Claerbout in 1970 for geophysical migration of seismic waves, in which the 
aim is to construct images of geological formations underground by downward 
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extrapolation of sound wave reflection data measured at the surface [7] For subsequent 
developments see [4,8,25] and recent volumes of Geophysics Second, they were intro- 
duced indepedently by Tappert in 1977 for underwater acoustics calculations, where 
the simulation of waves in the ocean can be speeded up greatly by assuming one-way 
lateral propagation away from the source [26] Further developments in this area are 
described in [22,24], the last of which contains many references, and in recent volumes of 
the Journal of the Acoustical Society of America Finally, one-way wave equations were 
proposed again independently by Engquist and Majda in 1977 [10,11], in extension of 
earlier ideas of Lindman [23], for application as approximate absorbing boundary 
conditions in numerical computations where artificial boundaries must be introduced to 
limit a computational domain For subsequent developments of this and related ideas, 
see [3,9,18] and recent volumes of the Journal of Computational Physics. 

Solving for f in (1 3), we get 

£ = ±tVI-t) 2 /t 2 , (15) 

where \T denotes the standard branch of the square root By (1 4), the plus and minus 
signs correspond to leftgoing and rightgoing waves, respectively 

f = + r Vl - y 2 /-! 2 v x < 0, (1 6a) 

f = - r Vl - T) 2 / r 2 <=3> v x > 0 (1 6b) 

Consequently the ideal one-way equation would be an equation with dispersion relation 
(l 6a), or equivalently, 


f = r Vl - s 2 , s — — = - sin0. 

r 


IDEAL OWWE 


(1 7) 


In the (£, y) plane, the curves of constant r determined by this equation are concentric 
semicircles Since Vl - s 2 is not rational, however, (1 7) is not the dispersion relation of 
a partial differential equation, but of an equation containing a pseudodifferential opera- 
tor Such an operator is not local in space or time, and cannot readily be implemented 
numerically, especially when the problem is generalized to allow variable coefficients. 

To construct practical one-way wave equations, therefore, we will approximate (1 7) 
with the aid of a function 


r(s) = 


Pm(«) 

?n (« ) ’ 
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where p m and q n are real polynomials of exact degrees m > 0 and n > 0 with no com- 
mon zeros We say that r is a rational function of exact type ( m,n ) The aim is for 
r(s) to be a good approximation to Vl - s 2 on [-1, 1], From any such approximation, 
we will derive a corresponding partial differential equation by replacing (1 7) by the 
dispersion relation 


C — rr(s), s — — — — sin0 

T 


APPROXIMATE OWWE 


( 18 ) 


For example, suppose r is the Taylor approximant r(s) = l-^-s 2 Substituting this 
for Vl - s 2 in (1 7) gives 




r C-T7>' 


or equivalently 


= t 2 - 



which is the dispersion relation of 


1 

= u tt ~ y«j ,y 


This is the original one-way wave equation, now classical, and is called the parabolic 
equation In the general case (1 8), we multiply both sides by q n (i] / t) to clear 

denominators, and get 


) = T<i Pm {— ), d = max{m, n +1} (19) 

T T 

This equation has the form 

P V, t) = 0, degree (P) = d, (1 10) 

where P is a homogeneous polynomial with real coefficients in three variables Thus it 
is the dispersion relation or symbol of a homogeneous partial differential equation of 
order d in x , y , t 

The most usual method for obtaining r is by Pade approximation [10,11] (Our 
example above is the Pade approximant of type (2, 0).) The disadvantage of Pade 
approximants is that although they are very accurate for s 0, they are inaccurate near 
the singularities at s = ±1, and as a result, the corresponding one-way wave equations 
behave poorly for 0 ~ y, Therefore it is tempting to consider alternative choices of 
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r Two reasonable candidates that we treat here are Chebyshev (L^) approximation 
and least-squares (L approximation (Certain least-squares approximations for one- 
way wave equations have been investigated previously [2,23,29] ) There are many other 
possibilities too, such as Chebyshev-Pade approximation [12], and we hope to compare 
the practical merits of various approximation schemes in a later work 

The purpose of this paper is to investigate the well-posedness of one-way wave 
equations Depending on the application, various well-posedness questions arise, of 
which we will consider two 

Initial value problem (IVP). A one-way wave equation could be applied as a 
partial differential equation in the domain t > 0, i,y£R, subject to initial data at 
t = 0 Well-posedness refers to the existence of a unique solution whose norm at t=t 0 
can be estimated in terms of the initial data. 

Initial boundary value problem (IBVP). In absorbing boundary condition 
applications, the domain is x, t > 0, y 6 R, and the one-way wave equation is applied 
as a boundary condition along x= 0 for (1 1) Well-posedness is now the existence of a 
unique solution whose norm at t = 1 0 and along x= 0 can be estimated in terms of the 
initial data 

Our mam result is that for most methods of approximation of Vl - s 2 , problems 
IVP and IBVP are well-posed if and only if (m,n) lies in one of three diagonals of the 
table of approximants, or two in the case of even approximants Our method consists of 
reducing each well-posedness question to algebraic criteria involving the function r , 
whose relationships to each other are studied systematically Here is an intuitive expla- 
nation of why the result comes out so simply even for the relatively complicated problem 
IBVP an absorbing boundary condition permits wave propagation in one direction, 
while lll-posedness amounts to the possibility of propagation in the other Thus absorp- 
tion and well-posedness are naturally related 

Analogous three-diagonal results have been obtained previously for other problems 
In the numerical solution of ordinary differential equations, the Ehle conjecture, proved 
by Wanner, Hairer, and N^rsett [30], asserts that a certain class of discrete approxima- 
tions is stable in precisely three diagonals More recently, Iserles and Strang have esta- 
blished a three-diagonal stability result for discrete approximations to hyperbolic partial 
differential equations [17] The connections between these results and ours are at present 
not understood 

Well-posedness results for one-way wave equations have been proved previously in 
[1,2,9,11,29,31], and elsewhere In particular, Engquist and Majda showed in [11] that 
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Pade approximation leads to a well-posed problem IBVP for all m= n and m= n+2, 
but an ill-posed one in case (4, 0) Thus another way to summarize our results is as fol- 
lows we show that the boundary conditions proposed by Engquist and Majda are the 
only Pade absorbing boundary conditions that are well-posed, and we extend an analo- 
gous conclusion to other classes of approximations and to initial value problems 

Two topics must be mentioned that are not discussed in this paper First, we do 
not consider the well-posedness of one-way wave equations as evolution equations in x 
rather than t , although this is their most common use in geophysics and underwater 
acoustics Second, nothing is said here about discrete approximations to one-way wave 
equations, which may be unstable even when the differential equation is well-posed We 
hope to repair these omissions in the future. 
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2. Statement of results and outline of the paper 

Our results can be summarized as follows Throughout, problems IVP and IBVP 
are based on the one-way wave equation derived by (1 8) from an arbitrary real rational 
function r of exact type (m, n) with r ( 0 ) 7^ 0 , 00 Also, the integer Xmn 1S defined by 

{ 0 if m+n is odd , 

1 if m+n is even ^ ^ 


r(s 1 

Theorem 1. Problem IVP is well-posed if and only if the zeros and poles of v ' 

s 

are real and simple and interlace along the real axis 

Theorem 2. Problem IBVP is well-posed if and only if the zeros and poles of lill 
are real and simple and interlace along the real axis, and furthermore r (s ) > 0 for 
* € [- 1 , 1 ] 

These results imply 

Corollary 1. Problems IVP and IBVP can be well-posed only if n < m < n +2 
Corollary 2. IBVP well-posed =^> IVP well-posed. 

Corollary 1 has the following converse - 

Theorem 3. If n < m < n+2, and if r(s) interpolates Vl - s 2 at m+n+l+Xmn 
points in (- 1 , 1 ), counted with multiplicity, then problems IVP and IBVP are well-posed 

The approximants one might consider in practice usually satisfy this interpolation condi- 
tion In particular, this is true of Pade, Chebyshev, and least-squares approximants 
(The least-squares approximants in question are required to be even, for as we will 
explain, the situation becomes more complicated otherwise ) This implies finally 

Theorem 4. The Pade, Chebyshev, and least-squares families of one-way wave 
equations yield well-posed problems IVP and IBVP in precisely two distinct diagonals of 
the tables of approximants, namely m = n and m = n+2 
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A remark should be made about odd values of m and n Usually, as in Theorem 
4, r is an even function, and the resulting one-way wave equation is symmetric about 
0 = 7T In this event m and n are even, so the case m = n+1 of Theorem 3 does not 
occur But there are also applications in which one wants a one-way wave equation 
tuned asymmetrically to be accurate near some angle 0 Q 7^ 7 r, and here, m or n may be 
odd We have therefore considered it important to derive results for arbitrary m and n 
rather than assume these numbers are even 

Here is an outline of the remainder of the paper Section 3 is devoted to investigat- 
ing certain purely algebraic properties of r that constitute the heart of our argument 
First, Lemma I establishes three equivalent statements to the interlace condition of 
Theorem 1, assuming r(0)> 0 Thus a fuller assertion than Theorem 1 would be that 
problem IVP is well-posed if and only if r or -r satisfies the conditions of Lemma 1 
Next, Lemma 2 establishes three equivalent statements to the condition of Theorem 2, 
problem IBVP is well-posed if and only if r satisfies the conditions of Lemma 2 One of 
these is an interpolation condition from which Theorem 3 will follow Finally, Lemma 3 
shows that a rational function r can interpolate Vl - s 2 in at most m+n +1+Xm» 
points in the plane This conclusion has been used already in the proof of Lemma 2, and 
will be applied again later to establish nondegeneracy of Pade and Chebyshev approxi- 
mants Because of it, Theorem 3 would be unchanged if it read at least m+n +1+Xmn 
points of interpolation 

Section 4 treats well-posedness of problem IVP by reducing it to one of the alge- 
braic conditions of Lemma 1 This proves Theorem 1 

Section 5 treats well-posedness of problem IBVP by reducing it to one of the condi- 
tions of Lemma 2 This proves Theorems 2 and 3. 

Finally, Section 6 shows that various families of even approximants to Vl - s 2 , 
including Pade, Chebyshev, and least-squares, satisfy the interpolation condition of 
Theorem 3 Together with Corollary 1, this proves Theorem 4 
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3. Lemmas on the rational function r 

Our first lemma identifies a class of functions that decrease monotonically along R 
except for simple poles Parts of this result have appeared elsewhere, for example in 
[1,2,31] 


Lemma 1. Let r be a real rational function of exact type ( m,n ) with 
r (0) 0,oo The following conditions are equivalent 

(a) If s is finite with Ims > 0, then r ^ is finite with Im r ( s ) < 0 

s s 

(b) r (0) > 0, and the zeros and poles of r ^ in C are real and simple and mter- 

s 

lace along the real axis 

(c) - can be written in the form 


s 

for some a , 6* , s k £ R, where < s 2 < * 


n+l 6 l 


t=i s ~ s k 

< s n+ i , 6 0 > 0, and b k > 0 for k > 1 


(3 1) 


Proof (a) =#> (b) We will show that if any of the four assertions of condition (b) 
fails, then (a) cannot hold First, if r (0) < 0, then r(s)/s behaves like a positive multi- 
ple of -1/s near s = 0, and necessarily maps te into the upper half plane for small 
enough e Second, if r(s)/s has a zero or pole «g6C-R, we can assume Imso > 0 by 
symmetry, and r(s)/s must take values near s 0 in both the upper and lower half- 
planes Third, if r(s)/s has a multiple zero or pole s 0 6 R, then r(s)/s must take 
values in both half-planes for s in a one-sided neighborhood of s 0 with lms>0 
Finally, suppose the zeros and poles are real and simple but do not interlace — there are 
two zeros s x , s 2 £R with no pole in-between, or the reverse. Then r(s)/s has a local 
maximum or minimum a at some point s* 6 (s x , s 2 ), so r(s)/s-a has a multiple zero at 
s* , and this is inconsistent with (a) as before. 

(b) =#> (c) Condition (b) implies that r(s)/s has exactly n+1 simple poles and 
n+2 simple zeros, hence order at most O(s) at s — oo Therefore its partial fraction 
decomposition has the form (3 1) for some a, b k , s k £ R with s x < < s„ +1 and 

with b k 0 for k > 1, and if k 0 is the index with s* 0 = 0, then b ko = r (0) > 0 Now if 

b k < 0 for some k > 1, there must be a pair b k , b k+l with b k b k+1 <0 But in this 
event there are an even number of zeros between s k and s t+1 , contradicting the inter- 
lace condition On the other hand if b k > 0 for all k > 1 but b 0 < 0, then r(s)/s has 
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constant sign and is therefore nonzero on (-00, Sj) and (s n+1 , 00) It follows that if the 
interlace condition holds, r(s)/s has only n real zeros out of n+2 zeros altogether, 
which implies the existence of a conjugate pair of zeros off the real axis, contradicting 

(b) 

(c) => fa). Trivial | 


The conditions of Lemma 1 are enough to ensure that r(s) does not interpolate 

-v/jT 

s 2 in the upper half-plane or at s = 0 For well-posedness of problem IBVP, this 
will have to be strengthened to exclude interpolation throughout [-1, 1] In what fol- 
lows, VW denotes the branch of the square root defined by branch cuts (-00, -1] and 
[1, 00) that takes the value 1 at s = 0. 


Lemma 2. Let r be a real rational function of exact type ( m,n ) with 
r (0) 7^ 0,oo The following conditions are equivalent: 

(d) r satisfies the conditions of Lemma 1, and moreover r(s ) > 0 for s E [-1, 1] 

(e) r (s) = -V 1 - s 2 has no solutions in C - (-00, -1) - (1, 00). 

(f) n < m < n+2, and r(s) = +Vl - s 2 has m +n +1+Xmn solutions in 
C - (-00, -1] - [1, 00), counted with multiplicity 

Proof (d) =^> (e). For Im s > 0, we have Im [-— \/l - s 2 ] > 0, while by condition 

s 

(a) of Lemma 1, Im r < 0 Therefore r(s) = -\/ 1 - s 2 has no solutions in the 

s 

upper half-plane, or by symmetry, in the lower half-plane. On [-1, 1] it cannot have any 
solutions either, since -V 1 - s 2 < 0 < r (s ) by assumption. 

(e) =^> (f) Squaring r (s ) = + Vl - s 2 yields the polynomial equation 

Pm{s)-(l~s 2 )q n \s) = 0, (3 2) 

which has exact degree 2 d= 2max{m, n +1}, hence has 2d solutions counted with mul- 
tiplicity, which can only lie in C - (-00, —1] - [I, 00) since (e) prohibits solutions at 
s = ±1 Each of these is a solution of r(s) = Vl - s 2 or of r (s ) = -V 1 - s 2 , but not 
both, so if (e) holds, there must be 2d solutions of the latter equation. Since 
2max{m,n+l} > m+n+l+x mn , with equality if and only if n< m< n+2, condition 
(f) now follows from Lemma 3, below 
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(f) => (d). Condition (f) can hold only if the estimates are sharp in the proof of 
Lemma 3, below, and this implies the interlace property (b) of Lemma 1; we omit the 
details. Moreover, r(s) must be positive on [-1,1], for otherwise r (s) = -Vl - s 2 
would have a solution there, and for n < m < n +2, this would put the total number of 
solutions of (3 2) above the limit 2max {m, n+1} | 

Finally, Lemma 3, whose proof we have just appealed to, limits the number of 
times r(s) can interpolate Vl - s 2 in the plane. The following argument is based on 
ideas suggested to us by Peter Borwein of Dalhousie University. See also Lemma 4 1 of 
[17], where Rolle’s theorem is used to limit the number of interpolation points on an 
interval Arieh Iserles has further indicated to us that an alternative proof of Lemma 3 
can be based on order stars [17,30] 

Lemma 3. Let r be a real rational Junction of exact type ( m,n ), and let L be the 
number of solutions of r (s) = Vl - s 2 m C - (-oo, -1) - (1, oo), counted with multipli- 
city Then 

L < m+n+l+Xmn- (3 3) 


Proof If there is a solution at s = ±1, its multiplicity is defined to be 1, since 
r(s) cannot match the infinite derivative there But we can dispose of this possibility 
by considering a new rational function r e (s ) = r (s + e) or r e (s ) = r (s +es ), depending 
on the signs of r(s) near the endpoints of [-1,1] For e sufficiently small and of 
appropriate sign, r £ = Vl - s 2 will also have L solutions, all of them in the open region 
C - (-oo, - 1 ] - [ 1 , oo). So without loss of generality, we can assume r (±1) 7 ^ 0 

The number L is the number of zeros of 

<f>{s ) = Pm (s ) “ Vl -s 2 q„ (s ) 

in C - (-00, - 1 ] - [ 1 , 00) For sufficiently large R , these must all lie in the region shown 
in Figure 1 Therefore <j) has winding number L along the boundary a-/3-'y-‘y-/3-a This 
implies that Im^(s) and Re<^(s) must each have at least 2 L zeros interlacing on this 
boundary In what follows we assume m > n+ 1 , and bound the number of zeros of 
Im<^(s), for m < n+ 1 , one proceeds analogously by bounding the number of zeros of 
Re <f)(s ) 
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Figure 1. Contour for counting zeros of r (s ) - \/l - s 2 

On /?, <f>(s ) behaves approximately like As m - tBs n+1 , where A and B are the lead- 
ing coefficients of p and q , respectively Therefore Im<^(s) has at most m + Xmn zeros 
on /?, where the Xmn results from the fact that if m and n are both even or both odd, 
the lower-order term iBs n+1 may bring about one extra crossing of the real axis On /? 
the count is the same On a, 7, 7, and a, we have Re <f>(s ) = p m (s ) and 
Im<^(s ) = -t \/l - s 2 q n (5 ), and therefore Im<£(s ) has at most 2n+2 zeros on these con- 
tours, including those at s = ±1 Adding these bounds together gives 

2L < 2m + 2x m „ + 2n + 2, 

which reduces to (3 3). | 
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4. Well-posedness of problem IVP 

Let r be a real rational function of exact type (m ,n), and consider the 
corresponding homogeneous partial differential equation of degree d— max{m,n+l} 
obtained from (1 8)-(1.10). Throughout this section, we assume r (0) =£ 0, oo, which is no 
restriction in practice since \/l - s 2 = 1 at s = 0 By (1 9), the condition r (0) 7^ 0 
amounts to the assumption that the coefficient of in P(£,ti,t) does not vanish, so 
that t — 0 is not a characteristic surface for the differential equation Here is the stan- 
dard result on well-posedness 

Well-posedness criterion. Problem IVP ts well-posed if and only if P (£, tj, t) = 0 
has no solutions with r] £ R and Im r < 0 


A partial differential equation that satisfies this condition is said to be hyperbolic 

For a precise definition of well-posedness and a derivation of this criterion, which is 
due originally to Garding, see Section 5 2(b) of [19] or Sections 12.3-12 5 of [15] The 
essential idea is Fourier analysis For suppose Problem IVP admits as a solution a mode 
(1 2) with f,i/6R and Im r < 0 By homogeneity, for any a > 0 there is another solu- 
tion with parameters (at;, at 7, ar) that grows at the rate e at l Imr l t and as a can be arbi- 
trarily large, the growth up to a fixed time t 0 cannot be bounded 

Here is the proof of Theorem 1 


r(s) 

Theorem 1. Problem IVP is well-posed if and only if the zeros and poles of — 

$ 

are real and simple and interlace along the real axis. 

Proof. Since the zeros of P (£,?/,•) are continuous functions of ij, m the well- 
posedness criterion it is enough to look for normal modes with y 7^ 0 By (1 9) and 
(1 10), since p m and q n have no common zeros, an equivalent well-posedness criterion in 
therefore that if £ £ R and y £ R - {0}, then 

t Pm ( n/r ) = 

V Vniv/r) V 


has no solutions with Im r < 0 -- or with Im r > 0, since the solutions come in conjugate 
pairs This is further equivalent to the statement that if v £ R for some 

s £ C - {0}, then s£R Contrapositively, if sGC-R, then r ^ £C-R Since 
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that v ' maps the upper half plane either into itself or into the lower half plane That 
8 

is, either r or -r satisfies condition (a) of Lemma 1 By that lemma, this is the same a 
saying that either r or -r satisfies condition (b), which is the assertion of Theorem 1 | 

As mentioned earlier, we have actually proved a stronger result: problem IVP is 
well-posed if and only if r or -r satisfies any of the conditions of Lemma 1 


5. Well-posedness of problem IBVP 

Again let r be a real rational function of exact type (m,n), and consider now the 
use of the corresponding one-way wave equation as a boundary condition at x = 0 for 
the wave equation (1 1) in the domain x > 0, y 6 R For simplicity we continue to 
assume r (0) 7^ 0, 00, as in the last section. . To determine well-posedness, we turn to the 
theory published by Kreiss in 1970 [20], specialized for the wave equation m the interior, 
and obtain: 


Well-posedness criterion. Problem IBVP is well-posed, if and only if 
P{£> T h T ) = Q an d ( 18 ) have no mutual solutions (£,77,7) 7^ (0, 0, 0) with rj E R, 
Im r < 0, and with £ belonging to the branch ( 1 . 6 b). 

Since r and £ are now complex, we have to be careful about what is meant by the 
branch (1.6b) For g £ R and Im r < 0, it means £ = -V t 2 - g 2 , the analytic function of 
r and g obtained by analytic continuation from the values £ = -r for g = 0 These 
values of £ satisfy Im £ > 0. For r 6 R, £ is defined by limits in the half-plane 
Im r < 0, and satisfies Im £ > 0 if |r|<|j/|,£6 [-1, 1] if | t | > | g | . 

For the definition of well-posedness and a derivation of this criterion, see [20] and 
[10,11]. Again, however, the essential idea is to look for ill-posed normal modes Sup- 
pose problem IBVP admits a solution (1 2) with Im £ > 0 and lmr<0 If Im £ > 0 
(decay as x— t-oo) but Im r < 0 (growth as f— ► 00), then by homogeneity, as in the last 
section, unbounded exponential growth as a function of ( is possible. If 
Im r = Im £ = 0 and (1 6b) holds, on the other hand, the solution is a plane wave with 
v x > 0 that radiates energy from the boundary into the interior, generating weaker but 
still unbounded growth See [13] for a presentation of the Kreiss theory from this 
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physical point of view, and [27] for the analogous discussion of stability of finite 
difference models. 

Here is the proof of Theorem 2. 


r(s 1 

Theorem 2. Problem IBVP is well-posed, if and only if the zeros and poles of v ' 

s 

are real and simple and interlace along the real axis , and furthermore r (s ) > 0 for 

s € [- 1 , 1 ] 

Proof If there is an normal mode (1 2) with r= 0, then f = ±trf, and by (1 9), 
the equation P (f, tj, r) = 0 reduces to 

± iAT d - Un 7 7 n+1 = BT d ~ m T) m , d = max {m ,n+l} 

for nonzero real constants A and B Whether or not m = n+1, this implies ij = 0 and 
thus £ = 0. Therefore we can assume r 7^ 0 Now for ^ G R and Im r < 0, the variable 
s — tj/t lies in C - (-00, 0) - (0, 00) Taking limits Im r— > 0 with r-b 0 amounts to let- 
ting s range over all of C, with points on the two sides of the cuts (-00, 0) and (0, 00) 
viewed as distinct By (1 6b) and (1 8)-(l 9), the well-posedness criterion now reduces to 
the condition that for s in this region, the equations f = -r Vl - s 2 and £ = rr(s) are 
never simultaneously satisfied. That is, the equation r (s ) = -V 1 - s 2 must have no 
solutions, where as m Section 3, the square root is defined by branch cuts (-00, -1] and 
[1, 00) and the value 1 at s =0 

On (-oo,-l) and (l,oo), r (s ) = -\/l - s 2 cannot be satisfied, because the left- 
hand side is real or infinite while the right-hand side is imaginary, finite, and nonzero 
Thus we have reduced the well-posedness crtierion to condition (e) of Lemma 2 By that 
lemma, an equivalent condition is (d), and by Lemma 1(b), this is equivalent to the 
assertion of Theorem 2 | 

Again, we have actually proved a stronger result problem IBVP is well-posed if and 
only if r satisfies any of the conditions of Lemma 2. In particular, by condition (e), it is 
well-posed if and only if r (s) interpolates \/l - s 2 at m+n+H-Xmn points in 
C - (-00, -1] - [1, 00) Therefore a sufficient condition for well-posedness is the existence 
of m+n+l+Xmrj points of interpolation in (-1, 1), and together with Corollary 2, this 
proves Theorem 3. 

If one or two points of interpolation lie at ±1 instead of in (-1, 1), problem IBVP is 
ill-posed This is a borderline case of weak ill-posedness, corresponding to a wave (1 2) 
that propagates tangentially to the boundary x — 0 
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6. Families of one-way wave equations 

In this final section we will examine various families of approximants r (s ) to 
Vl - s 2 on [-1, 1], which correspond to various families of one-way wave equations By 
Theorem 3, well-posedness is assured for n < m < n +2 if r (s ) interpolates Vl - s 2 at 
m+n+l+Xmn points of (-1, 1), and conversely, the approximants one might consider in 
practice will almost always satisfy this condition In fact, different methods of approxi- 
mation of Vl - s 2 can fruitfully be interpreted and compared as different strategies for 
the allocation of interpolation points in (-1, 1) -- hence of angles in (y, at which the 
ideal and approximate one-way dispersion relations coincide 

Until the end of the section, r is an even function First, we will show that 
s 2 can be interpolated in an arbitrary set of points symmetrically distributed in 
(-1, 1) Then we will consider Chebyshev, Pade, and least-squares approximation, and 
show that each leads to the appropriate number of interpolation points, thereby proving 
Theorem 4. (These methods of approximation are illustrative; we make no claim about 
what method may be best in practice.) In Chebyshev and Pade approximation, unique- 
ness implies that r (s ) is automatically even, for otherwise, r (-s ) would be a distinct 
and equally good approximation. In least-squares approximation approximation, unique- 
ness does not hold, so we have to impose evenness of r as an explicit constraint 


Symmetric interpolation in arbitrary points. The following construction has 
been pointed out to us by Philip Roe. For any even K > 2, let , • , ±s K / 2 be a 

set of K points in (-1, 1), counted with multiplicity, except that one of these pairs may 
be ±1 Set m = 2 and n == 2 [” -^-"j -2, so that m—n or m — n+2 and 

K = m+n+2 Let p be a nonzero polynomial of degree ~ that is zero at -y/l - s t 2 for 
each k , and set 


r (s) = 


p{t) + p(-t) 
{-p{t) + p{-t))/t’ 


( 6 . 1 ) 


where t = Vl - s " 1 Since the numerator is even as a function of t , it is a polynomial in 
s of degree m , and since the denominator is even as a function of t , it is a polynomial 
in 5 of degree n . Also, since | p (-t ) | > Ip (01 for t > 0 and p ' (0) 7^ 0, r (s ) can 
have no poles or zeros in [-1, 1] (we exclude the trivial case K = 2, s k = 1, r (s) = 0). 
Thus r (s ) = Vl - s 2 is equivalent to 


p (0 + pH) = -p (0 + pH )» 


that is, p (t ) = 0 In other words, (6 1) interpolates Vl - s 2 at the points ±s k , and by 
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Theorem 3, the corresponding one-way wave equation is well-posed 

An alternative and more transparent presentation of this construction is implicit in 
Section 9 of a paper by Higdon [14] If we write y/l - s k 2 = -cosd k , then p (f/r) = 0 is 
equivalent to 

if/2 

n (e + cos 0 k t) = 0, (6 2) 

*= 1 

which corresponds to the partial differential equation 

if/2 

II ( d x + cos 0 k d()u = 0, (6 3) 

*=i 

which is satisfied exactly by plane waves at the angles ±0 k . Conversely, given a one- 
way wave equation in the form (6 3), let us use (1 3) to eliminate all powers of £ higher 
than 1 in its dispersion relation (6 2) The result is a new dispersion relation that inter- 
polates Vl - s 2 in the points s k and is equivalent to a rational expression (1 8) For 
example, Higdon points out that the Engquist-Majda absorbing boundary conditions can 
be generated in this fashion from ( d x -d t ) K l 2 u =0, an observation that highlights 
their close relation to the boundary conditions of Bayliss and Turkel (see Section 2 of 

[ 3 ]) 


Now we turn to Chebyshev, Pade, and least-squares approximation In what fol- 
lows, let r denote a real rational function of type (M,N) but not necessarily exact type 
(M,N), that is, r has exact type ( m,n ) for some m < M , n < N Assume further, 
until further notice, that M, N, m, and n are even integers and that r is an even 
function The defect of r (as a function of type ( M, N )) is the integer 

8 — mm {M-m , N-n} > 0. (6 4) 

Also, let K be the number of zeros of Vl - s 2 - r(s) in (-1, 1) By definition K < L , 
where L is the integer of Lemma 3, and so by that lemma, since Xmn — 1 when m and 
n both even, we have 

K < m+n+2 (6 5) 


Chebyshev approximation. The Chebyshev approximant of type (M, N) to 
Vl - s 2 , denoted r^(s ), is defined by the condition 

II Vl-s 2 - rf 1N (s) 1 1 oo = minimum, 
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where || ||oo is the supremum norm on [-1, 1]. It is well known that exists, is 
unique, and is characterized by the following “equioscillation condition” [5,27] Let r be 
a rational function of type ( M,N ), exact type (m, n), and defect <5 Then r— rfj N if 
and only tf Vl - s 2 - r (s) attains alternating values ±||\/l-s 2 -r(s)|| 00 on some 
sequence of points -l < s 0 < < Sj < 1 with 

J > M+N+l-6 (6 6) 


Let r= rff N be the Chebyshev approximant to Vl - s 2 of type ( M,N ) Between 
any two equioscillation points there must be a zero, so we have K > J , and by (6 6), 

K > M+N+l-6 (6.7) 

Together with (6 5), this yields 

M+iV+1-6 < m+n+2, 

or by (6.4), 

M+N+l-6 < {M-6) + {N-8)+ 2. 

Therefore 8 < 1 Since 8 is even, however, we must have (5 = 0 Since K is even also, 
(6 5) and (6 7) now yield 


m—M, n = N, K = m+n+2 (6 8) 

(That is, the Chebyshev table for Vl - s 2 is normal except for the even-odd degeneracy 
it breaks into distinct 2X2 blocks of identical entries „ = r* +ln = r * n+1 
= r^ +1 „ +1 ) This establishes the hypotheses of Theorem 3, proving Theorem 4 for 
Chebyshev approximation 


Pade approximation. The Pade approximant of type ( M, N ) to Vl - s 2 , 
denoted rff N , is defined by the condition 

Vl~-s 2 - rh N {s) = 0(s m “ m ) as s — > 0 

Again one has existence and uniqueness (the former is trivial) For this particular case, 
rf fN is known explicitly and can be constructed from a continued fraction expansion 
[4,11,30], from hypergeometnc function identities [16], or from (6 1) Alternatively, to 
emphasize the analogy to Chebyshev approximation, one can reason by means of the fol- 
lowing equioscillation- type characterization [28]. Let r be a rational function of type 
(M, N), exact type ( m,n ), and defect 8 Then r = rjj^y if and only if 

Vl - s 2 - r (s) = 0(s J ) as s— ► 0 (6.9) 
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for some J satisfying (6 6) Since (6 9) asserts that r (s) interpolates Vl - s 2 at least J 
times at the origin, we again have K > J , and the argument leading to (6 8) goes 
through exactly as before (That is, the Pade table for Vl - s 2 is also normal except for 
even-odd degeneracy.) This proves Theorem 4 for Pade approximation. 

Least-squares approximation. Rational least-squares approximation, unlike its 
polynomial counterpart, has certain troublesome properties -- notably nonuniqueness and 
the possibility of local best approximations that are not global To begin with, let us 
drop the assumption that M, N, m , and n are even It is known that a least-squares 
approximant of type (M,N) to Vl - s 2 exists, so let us denote such a function by rj$f 
Here is the remarkable property that distinguishes this problem from Chebyshev and 
Pade approximation, proved by Cheney and Goldstein in 1967 [6]: rffy always has 6 = 0 
As a corollary, first pointed out by Lamprecht [21], cannot be even when M and N 
are odd, and so it cannot be unique We do not know whether it is unique when M and 
N are even. 

Nevertheless, certain conclusions about r$)i can be reached despite its lack of 
uniqueness. Since 5 = 0, one can show that \/l - s 2 - r must be orthogonal on [-1, 1] 
to all polynomials of degree M+N+l, and accordingly, must have K> M+N+l zeros 
in (-1, 1) (Cheney and Goldstein point this out as their final corollary ) Since 
K < m+n +1+Xmn by Lemma 3, we now reason as follows if M and N are both even 
or both odd, then m = M and n = N, but we are assured of only m+n + 1 interpola- 
tion points, not the m +n +2 required by Lemma 3. On the other hand if one of M and 
N is even and the other is odd, then the possibilities are m = M and n=N, 
m = M - 1 and n = N, or m = M and n = N- 1, and in each case, 
m+n +1+X mn = M+N+l, which guarantees a sufficient number of interpolation 
points 

We conclude that if one-way wave equations are derived by general least-squares 
approximation, problems IVP and IBVP are well-posed if M= N+ 1, but whether they 
are well-posed for M= N and M= N + 2 is uncertain 

Symmetric least-squares approximation. The situation returns to normal if, 
rather than permitting an arbitrary approximant of type (M , /V), we require r to be 
even. (Formula (6 1) will now be a computational help, for we can take the interpolation 
points as independent variables and vary them systematically to minimize 
|| v 1 - s 2 - r (s ) || 2 .) Assume again that M, N, m , n , and K are even The problem 
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is equivalent to that of approximating Vl - a on [0, 1] by a rational function of type 
(M/2, TV/ 2) in the least-squares sense with weight function <r -1 / 2 By an adaptation of 
the Cheney-Goldstein argument, we get 6 = 0 again, from which (6 8) follows. This 
completes the proof of Theorem 4 
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